— Original — 

Empirical Growth Curve Estimation Using Sigmoid Sub-functions 
that Adjust Seasonal Compensatory Growth for Male Body Weight 
of Thoroughbred Horses 

Tomoaki ONODA\ Ryuta YAMAMOTO^, Kyohei SAWAMURA^ Yoshinobu INOUE*, 
Akira MATSUI^ Takeshi MIYAKE^* and Nobuhiro HIRAli 

^Comparative Agricultural Sciences, Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, 
^The Japan Bloodhorse Breeders' Association, 4-5-4 Shinbashi, Minato-ku, Tokyo 105-0004, '^JRA Facilities 
Co. Ltd., 4-5-4 Shinbashi, Minato-ku, Tokyo 105-0004, ^Hidaka Training and Research Center, Japan Racing 
Association, 535-13 Nissha, Urakawa-cho, Hidaka, Hokkaido 057-0171, "^Equine Research Institute, Japan 
Racing Association, 321-4 Tokami-cho, Utsunomiya, Tochigi 320-0856, Japan 



Thoroughbred horses are seasonal mating animals, raised in northern regions or countries. J. Equine Sci. 

Foals born yearly in spring generally show a typical seasonal compensatory growth pattern. Vol. 22, No. 2 

in which their growth rate declines in the first winter and increases in the next spring. In this pp. 37-42, 201 1 

study, a new empirical adjustment approach is proposed to adjust for this compensatory 

growth when growth curve equations are estimated, by using 1,633 male body weights of 

Thoroughbreds as an illustrating example. Based on general Richards growth curve 

equation, a new growth curve equation was developed and fit to the weight-age data. The 

new growth curve equation had a sigmoid sub-function that can adjust the compensatory 

growth, combined with the Richards biological parameter responsible for the maturity of 

animals. The unknown parameters included in the equations were estimated by SAS 

NLMIXED procedure. The goodness of fit was examined by using Akaike 's Information 

Criterion (AIC). The AIC values decreased from 13,053 (general Richards equation) to 

12, 794 (the newly developed equation), indicating the better fit of the new equation to the 

weight-age data. The shape of the growth curve was improved during the period of 

compensatory growth. The proposed method is one of the useful approaches for adjusting 

seasonal compensatory growth in growth curve estimations for Thoroughbreds, and for their 

management during the compensatory period. Based on this approach, the optimal growth 

curve equations can be estimated also for female body weight of Thoroughbreds or other 

growth traits affected by seasonal compensatory growth. 
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Thoroughbred horses are seasonal breeders, and 
engage in their mating activities during early spring so 
that their foals bom 340 days later can take advantage 
of the milder temperatures and abundant forage in the 
next spring. Due to the seasonal mating, young 
Thoroughbred foals experience their first winter 
season almost simultaneously. They exhibit seasonal 
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compensatory growth that is characterized by a 
reduced growth rate in the first winter and an 
accelerated growth rate in the following spring. 

The seasonal compensatory growth is clearly found 
in the change of body weight or average daily gain of 
yearling Thoroughbreds especially in northern areas or 
countries [1 ] . Until now, there have been many studies 
investigating the seasonal change of horses' growth rate 
[2, 6, 13, 14] or determining the appropriate growth 
curve for the surveyed horses [8, 11]. However, 
empirical growth curve equations adjusted for the 
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Fig. 1. Scatterplot of the 1,633 male body weight (BWT) data of Thoroughbred coks (light gray dots). Thin black line 
indicates expected (actual) data averages of BWT, and the thick gray line indicates estimated Richards growth curve 
(i.e.. Equation 2). 



compensatory growth have not been estimated so far. 

Careful feeding management during the period of 
compensatory growth, based on useful standard growth 
curves, is important to secure the development of 
young Thoroughbred horses. The objective of this 
study was to estimate the single growth curve equation 
empirically adjusted for seasonal compensatory growth 
for male body weight of Thoroughbreds by using 
sigmoid sub-functions. In this study, the male body 
weight was only analyzed in order to illustrate an 
example for this approach. 



Materials and Methods 

A total of 1,633 male body weight (BWT) and age 
measurements, for 232 Thoroughbred colts, were 
collected by the Japan Bloodhorse Breeders' 
Association (JBBA) between 2005 and 2009. The data 
with extreme values over ± 4 residual standard 
deviations were preliminary analyzed by using Logistic 
growth curve equation, and removed from this study. 
The scatterplot of the weight-age data is shown with 
gray dots in Fig. 1. In the figure, the tendency of the 
seasonal compensatory growth can be recognized 
between about 300 and 600 days, which is the period of 
compensatory growth. 

Seven sigmoid growth curve equations (Logistic, 
Gompertz, von Bertalanffy, Brody, Richards, Bridges 



andjanoscheck; see [7] ) were fit to the weight-age data. 
The best-fit equation based on Akaike's Information 
Criterion (AIC) was the Richards equation [10] as 
follows: 



BWT 



{[.Q + Be-''') 



(Equation 1) 



In this "general" Richards equation (Equation 1), 
body weight {BWT, kg) is described as a function of age 
{t, day). The biological interpretations of the four 
parameters (A, B, k and m) have been discussed 
previously [13]. Briefly, A is the asymptotic limit for 
BWT as t approaches infinity; B is a scaling parameter 
that defines the degree of maturity; A; is a rate constant 
that determines the spread of the curve along the time 
axis; m is the point of inflection in the sigmoid curve in 
relation to age. 

These unknown parameters were estimated by the 
SAS NLMIXED procedure [12]. The estimated 
Richards equation using the weight-age data was as 
follows: 



BWT = 



575.0 



(1.0-0.945136-"'"""') 



1.0 (Equation 2) 



In Fig. 1, the thick gray line indicates the estimated 
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Fig. 2. Deviation of the expected (actual) data averages of the 
male body weight (BWT) from the estimated Richards 
growth curve (i.e., Equation 2). 



Fig. 3. Sigmoid curve with two different values of a. Black and 
gray lines are with a=10.0 and a=5.0, respectively. 



Richards growth curve Hsted on Equation 2, and the 
thin black Hne is a plot of the expected (actual) data 
averages of the BWT, where the averages were 
computed based on the monthly averages of the BWT 
data and their monthly daily gain. This thin black line 
is identical to that published in the Japanese Feeding 
Standard for Horses [4]. The deviation of the thin 
black line (expected data averages) from the thick gray 
line (Richards growth curve) is clear between about 300 
and 600 days. The deviation of the expected data 
averages from the Richards growth curve is shown in 
Fig. 2, where the largest deviations are recognized 
during the period of compensatory growth. 

To adjust for the deviation caused by compensatory 
growth shown in Fig. 2, a sigmoid sub-function was 
designed as follows. A general sigmoid function f(t) is 
expressed as: 



1.0 



(Equation 3) 



where e is the base of the natural logarithm, t is time, 
where -co < t<co, and 0 ^f(t) < 1. A sigmoid function 
is a monotonous increase function, and the shapes of 
the function depend on the a values. We chose a 
values of 5 and 10. The results around -1 < ^ < 1 are 
shown in Fig. 3. 

The subtraction between two sigmoid functions with 
different a values leads to a function having only a 
single wave (down and up wave) at ^=0. By using this 
subtraction and application of several modifications 




m 900 900 



Fig, 4. Constructed sigmoid sub-function f'(t) for adjustment of 
seasoned compensatory grovrth. 



around the t parameter, the following sigmoid sub- 
function /'(^j can be obtained: 



fit)= 



1.0 



1.0 



-10.0(?-432.0) 



1.0 + e 268.49 



-5.0(?-432.0) 
1.0 + e 268.49 



(Equation 4) 



and the shape of the sub-function is shown in Fig. 4. 
This sub-function curve crosses the x axis at 432 days 
with the wave length (i.e., distance between the starting 
and ending point of the wave) of 268.49*2. This sub- 
function has zero value when t = -co or oo. Only when t 
is around 432 days, do the non-zero /Y^) values appear. 
This characteristic of the sub-function is like the shape 
of the deviation shown in Fig. 2 especially for the 
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Fig. 5. Scatterplot of the 1,633 male body weight (BWT) data of Thoroughbred coks (light gray dots). Thin black line 
indicates expected (actual) data averages of BWT, and thick gray line indicates estimated Richcirds growth curve with 
the developed sigmoid sub-function (i.e., Equation 5). 



period of compensatory growth. Therefore, it was 
assumed that the f'(t) sub-function can empirically 
adjust the deviation caused by compensatory growth. 
In the original Richards equation (Equation 1) the 
maturity related parameter is B, and the newly 
developed f(t) sub-function can be used for the 
adjustment of the B parameter. 

Results and Discussion 

In general, the newly developed /'(O sigmoid sub- 
function can be added to the B parameter in Equation 1 . 
Specifically, this leads to the replacement of "-0.94513" 
by "-0.94513 + p*f'(t)" in Equation 2 where /? is a 
coefficient of f'(t). With this replacement and 
combination of Equation 2 and Equation 4, the growth 
curve equation acijusting for the seasonal compensatory 
growth was obtained as: 

BWT- ^ 

(l.0 + (-0.94513+0.3582/'(0)e"''°°'" ') ''''' 

r(A 1-0 1-0 

where J VV -10.0(r-432.0) -5.0(f-432.0) 

1.0 + e 268.49 l_0 + e 268.49 

(Equation 5) 

The optimal coefficient value of 0.3582 (i.e., the 



coefficient of f'(t)) and also the optimal wave length of 
268.49 were estimated again by using the weight-age 
data and SAS NLMIXED procedure. Equation 5 looks 
complicated but is a single function of age t in days. 
The insertion of f(t) in Equation 5 affects the BWT only 
when t is within about 432 ± 268.49 days. The growth 
curve of this new Richards equation (Equation 5) 
combined with the sigmoid sub-function is shown with 
thick gray line in Fig. 5, where the black line and dots 
are the same as in Fig. 1. As shown in Fig. 5, the shapes 
of the two lines were almost identical in the period of 
compensatory growth. By using this approach, the AIC 
value decreased from 13,053 (Equation 2) to 12,794 
(Equation 5), indicating the better fit of Equation 5 to 
the weight-age data than Equation 2. 

The deviation of the expected data averages from the 
new Richards growth curve adjusting for compensatory 
growth is shown in Fig. 6. The deviation in the period 
of compensatory growth is clearly reduced when 
compared to the case of non-adjusted (Fig. 2). These 
results suggest the usefulness of this proposed 
approach for handling compensatory growth typical in 
Thoroughbred horses. 

As Richards [10] noted, some researchers may 
consider that the parameter B is unimportant 
biologically. We can choose an alternative approach to 
adjust parameter k by the sub-function. If k was 
adjusted by the f'(t) sub-function, the combined 
equation became more complicated and difficult in 
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Age (days) 

Fig. 6. Deviation of the expected (actual) data averages of the 
male body weight (BWT) from the estimated Richards 
growth curve with the developed sub-function (i.e., 
Equation 5). 



numerical computation for parameter estimation 

because parameter k is the one used in the exponential. 
We aimed at simplicity for the combined equation, and 
chose B parameter to be adjusted in the f'(t) sub- 
function. 

Accounting for compensatory growth in 
Thoroughbreds, Staniar et al. [13] assumed that their 
growth rate consists of baseline and systematic 
deviation components. Even though they did not 
propose any unique growth curve equations handling 
compensatory growth, their idea shows a similarity to 
our approach of using both standard Richards function 
as the baseline and sigmoid sub-function as the 
systematic deviation. France et al. [5] showed a biphasic 
growth curve equation by combining two different 
equations with different phases (i.e., two different 
periods of age), and this may be an alternative 
approach to handle the compensatory growth because 
it has generally two different phases. Recently, Wan et 
al. [15] zndYmetal. [16] proposed their new sigmoidal 
growth curve functions but these equations seem not to 
be applicable for handling compensatory growth. 

Concerning developmental orthopaedic disease in 
horses, Donabedian et al. [3] investigated the 
relationship between nutrient intake and growth rate, 
and Raub [9] discussed the making of a durable 
equine. The standard growth curve equation handling 
compensatory growth would be useful for such 
experiments that relate secure development of 
Thoroughbred horses in winter and spring time. 

In conclusion, the proposed method in this study is 



one of the useful approaches for adjusting seasonal 
compensatory growth in growth curve estimations for 
Thoroughbreds. Based on this approach, the optimal 
growth curve equations can be estimated also for 
female BWT data of Thoroughbreds or another growth 
traits such as body height or entire width of chest that 
are considered to be also affected by seasonal 
compensatory growth. Further studies are also awaited 
for these traits. This approach is easily applicable to the 
general sigmoidal growth curve equation families such 
as those having biological parameters (e.g. A, B and k 
parameters in Richards equation). 
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